Lattice Induced Frequency Shifts in Sr Optical Lattice Clocks at the 10 Level 
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We present a comprehensive study of the frequency shifts associated with the lattice potential 
for a Sr lattice clock. By comparing two such clocks with a frequency stability reaching 5 x lO"'^^ 
after a one hour integration time, and varying the lattice depth up to Uo — 900 Er with Er being 
the recoil energy, we evaluate lattice related shifts with an unprecedented accuracy. We put the 
first experimental upper bound on the recently predicted frequency shift due to the magnetic dipole 
(Ml) and electric quadrupole (E2) interactions. This upper bound is significantly smaller than the 
theoretical upper limit. We also give a new upper limit on the effect of hyperpolarizability with an 
improvement by more than one order of magnitude. Finally, we report the first observation of the 
vector and tensor shifts in a lattice clock. Combining these measurements, we show that all known 
lattice related perturbation will not affect the clock accuracy down to the 10~^^ level, even for very 
deep lattices, up to Uo — 150 Er. 
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PACS numbers: 06. 30. Ft, 42.62.Fi, 37.10.Jk, 42.50.Nn 

Together with single trapped ion clocks, optical lat- 
tice clocks represent the future of frequency metrol- 
ogy. These two types of apparatus now outperform mi- 
crowave standards thanks to their five orders of magni- 
tude higher clock transition frequency and to a stringent 
confinement-based control of atomic motion. In trapped 
ion clocks, one takes advantage of the external charge 
to confine atoms in the relatively weak electromagnetic 
RF fields of a Paul trap [H. In contrast, optical lattice 
clocks [2] use neutral atoms and the trap based perturba- 
tions are orders of magnitude larger. Ion clocks currently 
exhibit the best accuracy at a level close to 10~^^ in 
fractional units [1, 0] • The best reported lattice clock ac- 
curacy is presently slightly above 10^^® with large room 
for improvement once one refines the control of the black- 
body radiation shift and of lattice related effects @] • In 
terms of frequency stability, lattice clocks are expected 
to outperform their single ion counterpart thanks to the 
larger number of atoms available in the experiment. 

The effects of the trapping potential are one of the 
main objects of study in lattice clocks. Indeed, at the 
theoretical minimal lattice depth Uq ^ 10 Er {Er be- 
ing the recoil energy associated with the absorption of 
a lattice photon) required for effectively cancelling mo- 
tional effects Q , the fractional energy shift of both clock 
transition levels is on the order of 10"^*^. In all prac- 
tical experiments to date, Uq is at least 5 to 10 times 
larger. Though the associated linear shift of the clock 
frequency is dramatically rejected by tuning the lattice 
to the magic wavelength Q , residual higher order effects 
have beenpredicted as potential limitations to the clock 
accuracy [8l4lOl|. In addition, polarization dependent ef- 
fects (vector and tensor shift, see below) have also been 
considered as possible serious issues for operating clocks 
with a 3D lattice with fermionic isotopes |lll4l3| . while 



this configuration is seen as the ultimate way to sup- 
press the frequency shift due to collisions between cold 
atoms Finally, a better control of the lattice induced 
effects would allow operating the clock at larger Uq than 
the theoretical minimum. This would make the clock 
operation easier, and more importantly, it would also al- 
low for a better control of several systematic effects. For 
instance, the shift due to cold collisions between atoms 
is expected to critically depend on Un via confinement 
strength and site-to-site tunneling 13l4l6l|. 

In this letter, we address these problems with a de- 
tailed study of all known lattice related effects. We show 
that they can be controlled at the 10"^'' level for a trap 
depth up to Uq — 150 Er. Based on these measurements, 
we also propose an optimal 3D lattice configuration for 
the fermionic *^Sr isotope. 

Three different lattice induced effects have been iden- 
tified. The first one is related to the non-scalar feature 
of the atom-lattice interaction, which results from the 
atomic hyperfine structure [1, 11 1 . This structure induces 
the minute dipole electric moment of the clock transi- 
tion allowing its laser excitation with extremely small 
linewidth. It also induces a small vector and tensor com- 
ponent in the atomic polarizability. In addition to the 
intensity dependence, this makes the light shift slightly 
dependent on the lattice polarization and geometry. The 
effect can be compensated for by tuning the lattice fre- 
quency, but this would be sensitive to possible long term 
variations of the laser's polarization and alignment as 
well as of the magnetic field direction. These effects have 
never been observed in a Sr lattice clock so far. The 
second effect, Hyperpolarizability, refers to a shift due 
to two-photon transitions. It was already pointed out 
in 2003 Q. The frequency shift scales as Uq and there- 
fore cannot be compensated for by a change of the lattice 
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frequency. A first experimental study sliowed that tiiis 
effect was small enough to not alter Sr clock accuracy at 
the 10-18 level for Uq = 10 Er ^. A better knowledge of 
hyperpolarizability is still required to ensure full perfor- 
mance at larger Uq . Finally, a subtle effect of higher order 
multipolar terms has recently been predicted [10[ . Dipole 
magnetic (Ml) and quadrupole electric (E2) interactions 
indeed lead to a shift of the internal atomic energy lev- 
els that is linear in lattice laser intensity. However, their 
spatial dependence in a lattice formed by the interference 
of several traveling waves does not match the one of the 
main dipole electric term. This in turn leads to a change 
of the oscillation frequency in the lattice potential wells 
and thereby of the spacing of the atomic motional states. 
This effect is expected to differ for the two clock states. 
The net M1/E2 frequency shift of the clock frequency 
scales as C/g^^ and again cannot be compensated for by 
tuning the lattice frequency. In Ref. fl(| it is argued that 
this effect could be as large as 10^^^ for Uq = 50 Er, a 
level at which the accuracy budget reported in ^ would 
have to be amended. 

Here, we measure the vector and tensor term with per- 
cent uncertainty. We improve the upper limit on the hy- 
perpolarizability effect by more than one order of mag- 
nitude, and finally perform the first experimental study 
of the M1/E2 effect showing that it is controllable to 
within IQ-i^ even for Uq as large as 700 Er- These mea- 
surements are performed by comparing two ^''Sr clocks 
(Sri and Sr2 in the following). Both clocks share the same 
clock interrogation laser at 698 nm 17, 3|. Fig.[T]shows a 
single scan over the resonance of Sri with a Rabi interro- 
gation of duration T = 250 ms. The resulting Fourier lim- 
ited linewidth of 3.2 Hz corresponds to a quality factor of 
Q = 1.3 X lO^"'; one of the highest ever obtained [l9l.[20|. 
The Allan deviation of the clock comparison is displayed 
in Fig.[T]and reaches a level of 7 x 10"^^ after one hour of 
averaging time. This record mid-term stability together 
with the possibility to operate the clock at extremely 
high Uo, up to Uo = 900 Er, are the key features of our 
apparatus for the study of lattice related effects. 

The lattice light at 813 nm is supplied by an all semi- 
conductor laser source. The setup consists of an extended 
cavity diode laser (ECDL), whose frequency is locked to 
a Fabry-Perot transfer cavity. The cavity itself is ref- 
erenced to the ^Sq — "^Pi atomic transition at 689 nm. 
The locked ECDL frequency exhibits long term fluctu- 
ations in the MHz range limited by the fluctuations of 
the cavity dispersion. Over a few hours, however, which 
is the typical duration of a measurement, the frequency 
is stable to within a few hundred kHz. The light from 
the ECDL is sent to two separate semi-conductor based 
amplifying systems providing up to 1.3 W to each of the 
lattice clocks. The high-power light is transmitted and 
spatially filtered through an optical fiber and coupled 
to a resonant cavity surrounding the atomic sample to 
create a ID lattice trap with depths up to 900£',-. Semi- 
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FIG. 1: Allan deviation for a comparison between the two Sr 
clocks. After an integration time of one hour, the t~^^^ trend 
reaches 7 x 10"^^, meaning that at least one of the clocks 
reaches a stability of 5 x 10"^^. Inset: single scan over the 
atomic resonance in Sri using a Rabi interrogation of 250 ms 
duration. The Fourier limited width is 3.2 Hz, corresponding 
to a quality factor of Q = 1.3 x 10^*. The contrast is close to 
90%. 



conductor amplifiers are known to produce non-negligible 
incoherent light spanning over tens of nm around the car- 
rier frequency, which is removed by an interference fil- 
ter of spectral width ^ 0.1 nm. We checked that the 
non-filtered incoherent light had a negligible effect on 
the clock frequency by changing the coherent/incoherent 
light intensity ratio for various lattice depths. 

We first consider the shifts resulting from the electric 
dipole (El) interaction. The clock frequency shift Av^^ 
can be split into scalar, vector and tensor contributions, 
which for a given \F, mp) state can be written (l8. [21| 



= (Ak" -I- Ak^tuf C4 • es + Ak'/3) Uq, (1) 



where /3 = (3|e ■ BbV - 1)[3to|, - F{F + 1)], 4 and cb 
are unitary vectors along the lattice wavevector and the 
quantization axis, respectivelty, e is the complex polar- 
ization vector of the lattice light, and ^efe ^ ie Ae* its 
degree of circularity. 

For each data point presented below, three to four se- 
quences with different lattice depths are interleaved to 
provide a determination of lattice related shifts. A fit of 
the motional sidebands provides the oscillation frequency 
in the trap as well as the longitudinal and transverse tem- 
peratures. From these, we determine the average trap 
depth Uq over the thermal distribution of the atoms in 
the trap. For each value of Uq, measurements with the 
symmetrical mp = ±9/2 states are alternated. Their 
degeneracy is lifted by applying a bias field Bq. 

The vector component in ([Ij is observed by measuring 
the frequency difference Zg = v(mp) — y{—mp) between 
two opposite mp states, which is insensitive to the scalar 
and tensor components. Its linear dependence on Uq re- 
flects the vector shift and its extrapolation to C/q = 
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FIG. 2: Linear component in Uo of the Zeeman splitting be- 
tween the mp = 9/2 and mp — —9/2 lines for different ori- 
entations of the quantization axis es. This orientation is set 
by adding a small magnetic field 6 parallel to the lattice, in 
addition to the usual orthogonal bias field Bq (Bq = 360 /^T). 
The scalar product Ck ■ Sb ~ ek ■ b/ Bq is deduced from the 
Zeeman splitting Za(Uo = 0) oc ^/fc^ -f B^. As expected from 
IH), the dependence is linear and its slope gives |9Ak''|. 



FIG. 3: Linear dependence of the average clock frequency 
on the trap depth. We vary the /3 parameter by changing 
the orientation of the quantization axis es, by probing the 
IttifI = 7/2 and \mF\ = 9/2 line and by switching between 
the two linear polarization eigen-modes of the lattice cavity. 
The lack of vector shift observed in Sri ensures that the po- 
larization of the lattice light is linear to a high degree. For /3 
= 0, we observe no dependence of the clock frequency on the 
trap depth, meaning that the lattice is tuned to the magic 
wavelength. 



gives the first order Zeeman splitting. As expected from 
the vector shift vanishes for a linear lattice polariza- 
tion, as well as for a bias field orthogonal to the lattice. 
We perform an accurate determination of Ak^ by oper- 
ating Sr2 with a circular lattice polarization (|^| = 1) for 
various bias field orientations, using Sri as a stable fre- 
quency reference. The results are displayed in Fig. [2] A 
linear fit to these data yields [23| 

IAk''] = {0.22±0M)Ez/Er. (2) 

Note that the vector shift does not directly affect the 
clock frequency, since it is rejected on average by alter- 
nately probing symmetrical Zeeman states. If the lattice 
polarization is linear, the vector shift also does not affect 
the bias field calibration, which must be done with per- 
cent accuracy to determine the frequency shift Ai/"? due 
to the second order Zeeman effect. We measured its the 
coefficient: 

^ -0.246(2) Hz/kHz% or = -23.5(2) MHz/T^ 

(3) 

The linear component of the dependence of the average 
clock frequency v = (^{mp) + v{—mp)) /2 on the trap 
depth Uo is (Ak* -|- Ak*/3)C/o, which gives access to the 
scalar and tensor polarizabilities. By varying the param- 
eter P on Sri while using Sr2 as a frequency reference 
(Fig. 131), we measure the tensor shift coefficient 

Ak* = (-0.0577 ± 0.0023) mHz/£;,.. (4) 

Its uncertainty is dominated by the knowledge of the 
atomic temperature, that determines the average power 
experienced by the atoms. This is the first observation 
of the tensor shift in a lattice clock. If the fluctuations 



of p are kept smaller than 1, the induced clock frequency 
fluctuations due to the tensor shift will be below 10~^^ 
for Uo — IQQEr- This would correspond to fluctuations 
of £*• cb at the 10"'^ level which can easily be achieved for 
any lattice polarization/bias field configuration. By op- 
erating the experiment in a configuration with a second 
order polarization dependence oi fi (|e • eep = or 1), /3 
can certainly be kept constant to within 0.1 or less. This 
confirms that the tensor shift will not limit the accuracy 
of Sr lattice clocks. 

Moreover, one can take advantage of the tensor shift to 
get rid of polarization inhomogeneities in a 3D lattices. 
This can be done by choosing different polarizations 
(hence different values of (5) for the lattice beams propa- 
gating in different directions and accordingly tuning their 
frequencies so that the total linear shift (Ak* -|- AK*/3)t/o 
is canceled everywhere in the lattice. The frequency dif- 
ference would wash out any polarization interference be- 
tween orthogonal directions, giving a perfectly homoge- 
neous effective polarization in the lattice. If for instance 
one uses je*- esp = i.e. (3 — —36 in two directions, and 
\e ■ eep = 1 i.e. /3 = 72 in the third one, the frequency 
difference is 300 MHz and the configuration is such that 
the /3 dependence on possible polarization fluctuations is 
second order for all the lattice beams. 

Making /? = gives an accurate determination of the 
magic wavelength of the lattice for which the scalar term 
Ak* in ((!]) is canceled. We find 

= 368 554 693(5) MHz, (5) 

which is in agreement with previous measurements [sj, 
with an accuracy improved by more than one order of 
magnitude. 
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FIG. 4: Left: hyperpolarizability extracted from the tensor 
shift data on Sri shown in Fig. [S] Right: hyperpolarizabihty 
obtained from a preliminary measurement of the vector shift 
on Sr2 with \^\ ~ 0.1 for various orientations of e*B. The latter 
exhibits a higher resolution due to a deeper trap. The grayed 
areas show the 1-a deviation from the weighted mean. 
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FIG. 5: M1/E2 dependence of the average clock frequency v 
of Sri. We observe no dependence of the clock frequency, but 
can put an upper bound on the coefficient. 



The hyperpolarizability is the Uq dependence of the 
average clock frequency v. Its coefficient was extracted 
by fitting the same data that gave the vector and tensor 
shift measurements with a parabola, and is plotted in 
Fig. m Though it is not resolved here, the hyperpolariz- 
ability shift is expected to depend on the lattice geome- 
try j22j . However, whatever this dependence, there exists 
at least one configuration where the effect is smaller than 
the weighted mean of the two datasets: 



(6) 



Its uncertainty is a factor of 60 lower than the previous 
evaluation reported in [9|], and is at the level of 10~^^ or 
lower if the lattice depth is kept at Uq < 150Er. 

For evaluating the frequency shift due to the M1/E2 
term, we use both the dependence on Uq and on the 
atomic motional state in the lattice. Since Ml and E2 
interactions modify the oscillation frequency in each lat- 
tice well, the shift of the clock transition for an atom in 



motional state \n) can be written 

A„^y^^ = C{n + 1/2)^0 



(7) 



Eq. [7] holds for a harmonic trapping potential, which is a 
good approximation for the large lattice depths and low 
atomic temperatures considered here. Experimentally n 
is varied by changing the atomic temperature T along the 
strong confinement direction while keeping the transverse 
temperature constant. The effective quantum number is 
then (n) = (exp(/iz/;/fcBr) — 1)^^, where vi is the oscila- 
tion frequency along the lattice axis. Measurements are 
performed by interleaving configurations with different 
lattice depth and temperatures, with Uq and (n) ranging 
from 50 to 400 Er and from 0.5 to 3.5 respectively. The 
data are displayed in Fig. \5\ giving the coefficient 

C = (0±0.31)mHz/V^. (8) 

The shift is not resolved here but gives a very stringent 
upper limit on the effect in optimal operating conditions. 
For the coolest configuration achived here ((n) — 0.5) 
and a trap depth of 150 Er, the uncertainty of the Ml /E2 
shift is smaller than 10^-'^^. We therefore conclude that 
this effect is much smaller than first feared and does not 
constitute a threat to the ultimate performance of Sr 
clocks. 

In conclusion, we have evaluated all relevant shifts due 
to the trap for a Sr lattice clock, showing that they will 
not thwart the ultimate performance of the clock. Sr 
clocks can even be operated at relatively large potential 
depths {Uq ^ 150 Er or more) while still keeping the lat- 
tice related perturbations below 10"^''. We also propose 
using the tensor shift for designing a 3D lattice configu- 
ration with no polarization gradients and a perfect can- 
celation of the linear shift. 
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